=3 



Study of Low Temperature Magnetic Properties of a Single Chain Magnet With 
Alternate Isotropic and Non-CoUinear Anisotropic Units 

Shaon Sahoo, Jean-Pascal Sutter and S. Ramasesha 
Abstract 

Here we study thermodynamic properties of an important class of single-chain magnets 
(SCMs), where alternate units are isotropic and anisotropic with anisotropy axes being non- 
coUinear. This class of SCMs shows slow relaxation at low temperatures which results from 
the interplay of two different relaxation mechanisms, namely dynamical and thermal. Here 
anisotropy is assumed to be large and negative, as a result, anisotropic units behave like 
canted spins at low temperatures; but even then simple Ising-type model does not capture 
the essential physics of the system due to quantum mechanical nature of the isotropic units. 
We here show how statistical behavior of this class of SCMs can be studied using a transfer 
matrix (TM) method. We also, for the first time, discuss in detail how weak inter-chain 
interactions can be treated by a TM method. The finite size effect is also discussed which 
becomes important for low temperature dynamics. At the end of this paper, we apply this 
technique to study a real helical chain magnet. 

Keywords Single chain magnets -Alternate isotropic and anisotropic units -Inter-chain interac- 
tion 



1 Introduction 

In the area of magnetism, the field of single chain magnets (SCMs) has attracted huge interest 
due to its potential use in dense data storage. Magnetism in one dimension is not only interesting 
from academic point of view, it has also created wide interest for application in technology. A 
huge body of research has already gone into understanding the physics of SCMs with the aim 
of designing systems with desired properties [[Tl|2l[3l|4l[5l. 

Though it is theoretically impossible to have a spontaneous magnetization at finite temper- 
ature in strictly one dimension, apparently there is no upper limit for the relaxation time of 
induced magnetization. In the field of SCMs, the focus therefore is to increase the relaxation 
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time by combining two modes of relaxations, namely, dynamic relaxation (creation and move- 
ment of domain wall(s) along chain) and thermal relaxation (hopping across anisotropy barrier). 
This is expected to be achieved by synthesizing a spin chain with anisotropic units. While syn- 
thetic chemists are exploring different chemical possibilities, theoreticians are developing new 
techniques to study these compounds. The theoretical studies in return help to understand the 
paradigms involved in increasing the relaxation time. 

Modeling a SCM has to take into account the chemical and geometrical structure of chain 
units that the system possesses |3. With quantum units, exact numerical diagonalization 
technique is not feasible, due to large size of the system, particularly so at low temperatures as 
correlation length grows exponentially as the temperature is lowered and the system size to be 
studied should at least have dimension of the correlation length. While a large number of SCMs 
have been synthesized over the years, here we will study an important class of SCMs in detail. 
In this class of SCMs, units are alternately isotropic and anisotropic. Furthermore, anisotropy 
axes are not coUinear and since anisotropy is large and negative, the anisotropic units can be 
thought of as canted Ising spins at low temperatures. It is worth noting that, in this situation had 
the anisotropy axes been coUinear, we could have used a simple one dimensional Ising model as 
only one component of the spins of the isotropic sites would be coupled to the anisotropic units. 
But the class of SCMs we wish to model, has non-coUinearity of the anisotropy axes, therefore 
different components of the spin of the isotropic sites would be involved at different sites; as 
a result, we have to treat the isotropic spins quantum mechanically. Here we note that though 
all the interactions are quantum mechanical in nature, some of them become insignificant at 
low temperatures due to highly preferred spin-direction of some sites. Herein, we model these 
type of systems exactly by transfer matrix method and study their low temperature statistical 
properties. We also discuss the issue of interchain interactions explicitly, which has been a long 
pending issue in this field O. 

There are some related theoretical/numerical studies which deserve a brief mention here. 
Alternative classical-quantum (or anisotropic-isotropic) spin systems have been studied previ- 
ously [l8l|9l[T0l[IIl[l2l[T3l[T4l. There, classical nature of the spins is assumed from correspon- 
dence principle, due to large value of the spin S, of the quantum spins. In an another work, a 
helical spin chain with alternately classical spin (Ising-type) and quantum spin (coming from 
organic radical) spins is studied [|T5l . But here the quantum spins are reduced to classical spins 
by approximating Ising direction by averaging directions of their neighboring classical spins. 

Before we present our method in Sec. 3, we present in the next section (Sec. 2) a brief 
discussion of a real system which can be modeled to be in the class of SCMs we are interested 
here. This real example should clarify the present modeling issues and provide the motivation 
for the model studied in this paper. In Sec. 4, we apply our technique to the real system and 
present our results. We conclude our work in Sec. 5. 
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Figure 1 : One structural unit of Fe2Nb chain is shown here. CO is the angle made by anisotropy axis of 
in-chain Fe-ions with chain-axis. 5 is the angle made by anisotropy axis of off-chain Fe ion with z-axis 
of associated A'^Z^-ion. a is the angle made by local xy-plane of Nb with the plane formed by all alternate 
in-chain Fe-ions (here Fe\ and Fej) along with their anisotropy axes. Alternate classical and quantum 
units are also indicated here by dotted boxes (see Sec. 3). 

2 Description of a Real System 

Recently a compound, we refer to as Fe2Nb, whose chemical formula is [{{H2O) Fe{L)} 
{Nb{CN)g} {Fe{L)}]^ has been synthesized lfT6l . Its crystal structure studies reveal that it 
consists of topologically quasi-one dimensional spin chains with helical geometry. This system 
shows SCM character of slow relaxation and large DC susceptibility at low temperatures. The 
back bone of each chain consists of alternate Fe and Nb ions with each Nb ion connected to an 
additional off-chain Fe ion (see Fig[il). Nb ions have spin 1/2 while Fe ions have spin of 2. The 
coordination around Fe{U) is the unusual hepta coordination which leads to large anisotropy 
of the Fe spins. All the intrachain interactions are antiferromagnetic and much stronger than 
interchain interaction. Since the Fe ions have large and negative (easy-axis) anisotropy, we can 
assume the Fe spins to be Ising-type, with non-coUinear spin orientations due to helicity of the 
chain. As the spin axes are non-coUinear, different components of Nb spins are involved in in- 
teractions at different sites along the helix and therefore must be treated quantum mechanically. 
A careful examination of the chain structure, shows that two chemical units form a geometri- 
cal unit which are used as basic blocks in our transfer matrix method. To define this complex 
structural unit, we introduce three structural parameters, CO, a and 5, which are defined in Fig 
(dl). We note that, here a quantum unit is formed by two A'^^ ions connected by one in-chain Fe 
ion and each of the Nb ions being connected to an off-chain Fe ion. These quantum units are 
connected by in-chain Fe ions, which behave like classical spins. Therefore, this system falls 
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Figure 2: Schematic diagram of a chain with alternative quantum and classical units 



into the class of SCMs we discussed in the previous section. In Sec 4, we first discuss in detail 
the application of our technique to solve the spin model for the SCM system. This is followed 
by presentation of results and their comparison with experimental studies. 



3 Description of the Technique 

Let us consider our system formed by identical quantum units connected by classical (Ising- 
type) spins (see Fig [21). Let 'H.f^{Qi) be the Hamiltonian associated with i-th quantum unit 
with the set of quantum and classical spin operators Q, . Let the Hamiltonian for the interaction 
of the ?'-th quantum unit with classical spin, S,-, to the left, be H^'^(qf', S,) with qf the quantum 
spin operator of the i-th quantum unit involved in interaction to the left. We can similarly define 
H-' (qf , S/+i) as the Hamiltonian for the interaction of the i-th quantum unit with classical spin, 
to the right. In the presence of (uniform) external magnetic field B, there will be Zeeman 
interaction term —g^B^z,iB, where Sz,/ be the Z-component of spin for the z-th block, and this 
term is here interpreted as sum of the projection of spin for the i-th quantum unit and half of 
that for both i-th and (/ + l)-th classical spins along the global (laboratory) Z-axis. Here Sz./ is 
defined with weight 1/2 for the classical spins to avoid double counting. We now can write our 
total Hamiltonian of a chain as, H = YJi=\ H/, where, 

= Hf'^(qf,S,)+HP^(Q,0+H;'^(qf,S,-+i)-gA'sSz,,5, (1) 

In the last term, the g-factor is assumed to be an average over all the corresponding g-factors 
of different spin-sites and /jg is the Bohr magneton. We note that, the Zeeman interaction term 
has three different parts: 5Hf'"(S/,5), Hf'"(Q/,5) and 5H^"j(S,+i,5) corresponding to the one 
associated with z-th classical spin, z-th quantum unit and (z + l)-th classical spin respectively. 
To avoid double counting, as before, Zeeman interaction part for z-th and (z+ l)-th classical 
spins have been taken with weights 1/2. Partition function for the chain (having N quantum and 
classical units) at temperature T can be written as, 

<2yv(P,5)= Tre-PH (2) 

where, 'Tr' means trace and (3=1/ {ksT), kg being the Boltzmann constant. Now let us consider 
direct product basis of the whole chain. This can simply be written as, | ■ ■ ■ , o,, Qi,- ■■), where 
|a,) and \Qi) represents convenient basis for z'-th classical spin and quantum unit, respectively. 
Eq ^ can now be written as, 

Qn{^.B)^ £ (■■■,o„avk"'^"|---,cr,-,a,---) (3) 
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Here the sums run over all possible configurations {a, Q} of all the classical spins and the 
quantum units. Now since, and commute, we can rewrite Eq (|3]) in the following form: 




QNi^.B) = £ (■■■,o„avin^"'^"'l---'<'^-'av) 

••■,a,v) (4) 

Here, Qj is the set of configurations for i-th quantum unit. Note, the quantity in parentheses, 
52(2/k '^^' iQi)^ is apurely classical operator (containing S/ and S^+i only) as we have summed 

Qi 

over quantum variables of the unit. Denoting it as T, enables us to rewrite Eq Q as, 

N 

Qn{^.B) = £(...,a,vinTH---,o.v) (5) 

{a} 1=1 

This form is familiar to us, with T, being the transfer operator. Introducing the identity ^ |o/) (o, | 

between successive Ts, will enable us to write the Partition function as a trace of A'^-th power of 
a transfer matrix, P, 

Qn{%B) = trP^ (6) 
where, elements of the transfer matrix P are given by, 

= {<S^W<S^+l) (7) 

Since we have taken classical spins to be Ising-type (which can have only two values, ±1), 
the P matrix will be a 2 x 2 matrix. Diagonalizing this matrix is trivial; let X+ and X- be 
the two eigenvalues of the matrix, then they can be expressed in terms of four elements of P 



as: A, 



± 



{pn+ Pii) ± V (j^ii -7^22)^ + 47712^21 



. In terms of eigenvalues, the partition 



function can now be written as, 

Qn{^,B) = Xl + X^l (8) 

Note both the eigenvalues are function of P and B. For thermodynamically large A'^, one can only 
take the larger eigenvalue (here A,-|-), but for finite A'', one should take both of them to evaluate 
the partition function. %r as a function of T (% being susceptibility), can be obtained as, 

where, A^^ is the Avogadro's number. In Eq Q, 2nd order differentiation can be done easily, 
one can even get simple closed analytical form for small applied magnetic field B. In general, 
one can do simple numerical differentiation to get xT value as a function of temperature T. 
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We now need to discuss a little more about T/. If we take basis \Qi) to be eigenstates of 
H,-, then T,- will just be sum of exponential of those eigenvalues. Let by diagonalizing it (if 
quantum unit has unconnected spin- 1/2 quantum spins, we can do it trivially and get analytical 
expression) we get eigenvalues L^ , where j runs from 1 to J (number of basis vectors or dimen- 
sionality of configurational space for the quantum unit). Note these eigenvalues are functions 
of only S/ and S,+i, the classical spin operators. So we can write T,- in the following way: 

= £(a-k-P«'ia) 

Qi 

= ie-^^J (10) 

7=1 

In our case (Fe2Nb), quantum unit has five spins (3 Ising spins and 2 quantum spin-l/2s), hence 
d = 32 and T,- will involve summation over the 32 eigenvalues. 

3.1 Consideration of Interchain Interaction 

Interchain interaction is an important issue in this field of SCMs; in real systems it is not possible 
to separate chains to eliminate this interaction. With increase in this interaction strength, a 
one dimensional system like SCM gradually transforms into a higher dimensional system with 
remarkably different properties. Even a weak interchain interaction can change properties of a 
chain -its presence can be seen both in static (like xTvs.T plots, hysteresis curves) and dynamic 
(like, Cole-Cole plots) measurements. Unfortunately, there is no detailed and explicit discussion 
on how one can treat this interaction theoretically. In fact, this issue has been highlighted in a 
paper by Miyasaka et al. [O, where it is posed as an open problem. In this section we discuss 
interchain interaction in detail. 

In low-dimension, the quantum fluctuations are generally large. However SCMs are not 
strictly one-dimensional, and we have weak interchain interaction. It is customary to deal with 
interchain interactions within a mean field approximation. Here a widely used formula for 
corrected (or modified) susceptibility due to the interaction is Xm = Xo/(l + zJ'Xo/NAg^ijjj), 
where Xo is susceptibility of the chain without interchain interactions, z is the number of nearest 
neighbors and f is interchain interaction strength (zf together is called mean field parameter) 
ifTTl [m [m. However we wish to emphasize that, this formula is for vanishing applied mag- 
netic field. So in case of finite magnetic field and at very low temperatures this formula is not 
appropriate. In addition, for a chain, unlike for molecular systems, Sz can be very large; so it 
is not clear whether it is valid to make series expansion of exponential function of the Zeeman 
and interchain interactions and retain terms only up to first order. Here we discuss a mean- field 
approach where the interaction is dealt by transfer matrix technique and can be used to study 
the case when magnetic field is finite. We also wish to make it clear that, there has to be a 
factor A'' (number of chemical units in a chain) multiplying zf in the above formula of modified 
susceptibility for vanishing magnetic field. To the best of our knowledge, this is the first attempt 
to deal with interchain interactions by the transfer matrix technique. 
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If the interchain interaction is weak, we can bring in the effect of environment by placing a 
chain in the mean field generated by neighboring chains. Let us now consider {Sz) be the expec- 
tation value of Z-component of spin for each chain and zJ' be the interchain interaction (mean 
field) parameter, then the interaction energy can be written as: zJ'Sz{Sz), with Sz being the 

A' 

Z-component of total spin for a chain. Now since Sz = ^ Sz,; with Sz,/ being the Z-component 

;=1 

N 

of spin for z-th block, we can rewrite the interaction energy as '^{zJ'Sz.i{Sz)), which can be 

/=i 

viewed as sum of interaction energy of each block of a chain with a mean field generated by 
neighboring chains. Now if we add this mean field block interaction energy term, zJ'Sz,i{Sz), 
to H,- of Eq ([B, we get, 

= H;'^(qf,S0+HP^(Q,-)+Hf^(qf,S,-+i)-^A'fiSz,/5 + z/Sz,/(5z) 

= H;'^(qf,S,) +hP^(Q,0 +H;'^(qf ,S,-+i) -^A'sSz,/5' (11) 

where, B' = B - ^{Sz). We can then follow the usual transfer matrix technique and get a 
partition function for the chain in the presence of a modified magnetic field B' . Since, (Sz) = 
_l_9in_a^Ml^ we see that, {Sz) has to be solved self-consistently by an iterative procedure. 
We can use the following self-consistent scheme developed by first noting that, zJ' being very 
small, internal field, AB = — — {Sz), is a slowly varying function of external applied field B 

resulting in | ^ | ^ 1 or 1 1|- 1 ^ 1 . Using the conditions and making the Taylor series expansion 
of In Qn{^,B') about B, we get, 

ain 2^(13,50 



gA'Bp {Sz) 



dB' 

din Qn{%B) f zJ'\dhnQN{^,B) 



dB \giuBj BB^ 



{Sz) + 



1 f zJ'Yd^\nQ!,{^,B) 



9\ y ;)R3 (5z)'---- (12) 
2\gHBj oB^ 

Iteration can be started with {Sz) = ^^g^'^^ , without interchain interaction, finally lead- 

ing to a converged {Sz) value. The self-consistent result correct to second order in zJ' is given 

by, 

s,b{Sz) - m(M)-(^)m(M)^ + 



zJ' 

2 2 
§ n 



dB 

- 1 ^d^M({\R)' 

(13) 



where, M((3,5) = jg the magnetization in the absence of interchain interaction. 

Going to higher order in zJ' is possible but generally not necessary. Obtaining modified sus- 
ceptibility (%m) from Eq (fT3l) is straight forward. In the vanishing field regime, we can get 
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the widely used fomiula for by retaining terms up to first order in zJ' of the Eq (fT2l) and 
subsequently solving for {Sz). This gives, 

1 9ln 

Now in this regime, modified susceptibility is %m = ^gPB{Sz) and susceptibility without in- 
terchain interaction is %o = ^g^^^^ • "'"'^^^ directly leads to: 

Xm = - Sly (15) 

Note the difference of factor of A'^ with this formula and the one generally used in this field 
[fT7l[T8l[T9l . Appearance of this factor in Eq (fT5l) is easy to understand: while susceptibility is 
measured for Avogadro number (Na) of chemical units, in the transfer matrix method we com- 
pute partition function for (chemical) units; therefore, we need to normalize the susceptibility 
computed from the partition function. The factor NzJ' can be given a physical interpretation. 
Let us assume that a unit in a given chain interacts only with the corresponding units of the z 
neighboring chains and p be the total effective coupling strength, then the associated interaction 
term will be pSz,;%^, where is the expectation value of Z-component of spin for a unit. 
The total mean field interaction term is therefore given by L,\pSz,; -^. Comparing this with 
conventional chain-chain total interaction term Y.iZJ'Szj{Sz), we get p = NzJ'. Thus, while zJ' 
is the coupling strength of a chain with its z neighboring chains, NzJ' is the coupling strength of 
a given unit with the corresponding units of the z neighboring chains. We will report the latter 
parameter (Nzf) for our study of real system. 

This way of introducing interchain interaction in computing magnetic susceptibility is fairly 
general and can be applied to any class of SCMs. Here we have presented a formula for modified 
susceptibility which is very similar to the one used for molecular systems and generally adopted 
for SCMs. The higher order self-consistent method (Eq[T3l) is more general and applicable even 
when the magnetic field is not vanishingly small. 



4 Application to the Fe2Nb system 

In this section we will apply our technique to the system Fe2Nb which is already introduced 
in Sec. 2. All the Fe — Nb couplings are supposed to be antiferromagnetic and are assumed 
to be of same strength (7). We assume that the strength of axial anisotropy constant D to be 
the same for each Fe ions. As we already mentioned, since this D is known to be large and 
negative, we treat the spins on Fe ions as two state Ising spins along their respective anisotropy 
axis. We also further assume the g-factor to be isotropic and has same value for all the ions. 
Three parameters are defined in Fig ([T]), namely, (O, a and 5, to describe complex structural unit 
of the chain. In Eq ^ we need to use average partition function, averaged over all possible 



8 



orientations of a chain to compare with experiments done on a powder sample. The averaging 
can be done by noting that a chain orientation in space can be defined by three angles, namely, 
0, (|) and \|/. Here and (|) are usual spherical polar coordinates to define the orientation of the 
chain axis and \\r is measure of rotation of the chain about its own axis. Our partition function 
will depend on (a) temperature, (b) magnetic field and (c) all three structural parameters (co, a 
and 5)for a single chain, it will also depend on orientational angles (0, (|) and \|/) with respect to 
laboratory frame. The partition function can be averaged over all orientations of the chain axis 
to obtain an average partition function 2Af((3,5), given by, 

pTi p2iz p2ii 

Qn{%B)= / / QN{^,B;Q,(^,\\f)^mQdQd(^d\\f (16) 
je=o>/(t)=o>/\|/=o 

where, Qn{^,B; 0, (|), \|/) is the partition function of the chain for a particular orientation (0, (|), \\f). 
We have ignored a normalization constant, as it will not affect our result. Integrations in Eq (fT6l) 
can be done by simple numerical methods (analytically it may not always be possible). So in 
this particular problem, Eq ^ must be replaced by. 

Note, we have also replaced N by 2N, as one structural unit consists of two chemical units and 
susceptibility from experiment is quoted for an Avogadro number of chemical units. 

The Hamiltonian for z-th quantum unit consisting of Nbi , Fe[ , Fe^-, Nb2 and Fe'2 site spins 
and connected to neighboring quantum units by Fei and Fe2 site spins respectively, can be 
written as: 

-^g/jsS^'^^.B-^glJBS^'^-B 



;S^^i ■ B + S^'^'i ■ B + ■ B + S^*2 . B + . ^) (18) 



Here S is a vector spin operator for the specified ion. All these spin operators are defined relative 
to their corresponding local coordinate systems (see Fig[il). For example, z-axis of a Fe ion is 
taken along its easy anisotropy axis. First and second terms in Eq (fTSi) are interaction terms 
between quantum unit and its neighboring classical units. The 3rd term within brackets includes 
all the interactions within the quantum unit. The 4th, 5th and 6th terms are the Zeeman terms 
associated with classical spins with appropriate weight and the quantum unit respectively. All 
these terms can be expressed explicitly with the help of structural angles or parameters (o), a, 5) 
and polar coordinates of the chain (0, (|), \\r). Table gives the expression for each of the terms 
in the laboratory frame assuming that the magnetic field is applied along Z-axis of the laboratory 
frame. Knowing all these interaction terms, now we can form the transfer operator T for the 
z-th unit as, 

T,- = £(a-|e-PH,'|2^.) 

{Qi} 

{Qi} 
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Table 1: Expression for all 13 interaction terms appearing in Eq ( fTSl ) is given below. Treating the Fe 
spins as Ising spins with anisotropy along z direction, we replace the spin by 20^*^', where a^*^' 
is the Ising- variable with only possible values ±1, the factor 2 is to account for the spin of the 
Fe ion which is 2. 



Term 


Expression 




2o^^i (Sx cosco cosa + sinco + sf''' coso) sina) 


gFe2 . '^Nb2 


2<5 ~[Sx coscocosa — Sv sinro + Sj; - coso) sina) 


^Fe\ . ^Nb, 


2o^^i(Sf' cosS + Sf' sin5) 




20^*^' (S^^' cosco cosa — S^^' sinco + sf''' cosco sina) 


gFfic . '^Nbi 


2(5^^'= {Sx'^- COSO) cosa + Sy''- sinco -1- ''^ cosco sina) 


^Fe'2 . ^Nb2 


20^4 (Sf 2 cos5 + Sj"^ sin5) 


S^^i .B 


20^*^1 ( cosco COS0 + sinco sin\|/ sin0)5 


^Fe2.§ 


20^*^2 cosco COS0+ sinco sin\|/ sin0)5 


^Nbi 


[Sx ' (cosa COS0 — sina sin0 cos\|/) + Sy ' sin\|/ sin0 
+sf (sina COS0 + cosa sin0 cos\|/)]5 


S^^'i ■ B 


20^*^1 (sin5 sin\|/ sin0 + cos5 sina cos0 + cos5 cosa cos\|/ sin0)5 


S^*^' ■ B 


20^*^' (cosco COS0 — sino) sin\|/ sin0)5 


^Nb2.g 


[Sx''^(cosa COS0 + sina sin0 cos\|/) — S^*^sin\|; sin0 
+sf (sina COS0 — cosa sin0 cos\|/)]5 


S^'^z . B 


20^*^2 (—sin5 sin\|/ sin0 + cos5 sina cos0 — cos5 cosa cos\|/ sin0)5 
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Expressions for a, b, c, d, e, f and h, which are functions of classical spins involved in the 
problem, are given in Table Here \Qi) is the direct product basis of the spins of Nbi, 

Table 2: Expression for variables appeared in Eq ( [T9l ) 



Variable 


Expression 


a 


— (3[27(a^^i + a^^'')cos(o cosa — g/j55(cosa cosG — sina sin0 cos\|/)] 


b 


-|3[2/{(o^^i -a^'^Osinro + a^^'isinS} -gA's^ sin\|/ sinB] 


c 


— |3[27{(o^^i +a^'^'')cos(0 sina + o^^'icos5} — g/JBjB(sina cosG + cosa sinG cos\|/)] 


d 


— P[27(a^^'' + o^^2)cos(o cosa — g/js5(cosa cosG + sina sinG cos\|/)] 


e 


-p [27{ (a^^-^^ -a^^2)sin(o + a^4sin5} +gA's5 sin\|/ sinG] 


f 


-|3[27{(o^^'^ +o^^2)cos(0 sina + o^^'2cos5} -g/Js5(sina cosG - cosa sinG cos\^)] 


h 


|3g/js5[(o''^' + a^'^-){cos(}i cosG + sino) sinG sin\|/) 

+2o^^i (sin5 sinG sm\\r + cos6 sina cosG + cos5 cosa sinG cos\|/) 

+2o^^'^(cos(0 cosG — sinco sinG sin\|/) 

+2a^^2(— sin5 sinG sin\|/ + cos5 sina cosG — cos5 cosa sinG cos\|/)] 



Fe[, Fee, Nb2, and Fe'2. If we take the basis for Nb\ to be eigenstates of (aS.^ ' -\-bSy ' + 
cSf*') (the quantum mechanical operator in the exponent of the first exponential in Eq (fT9l)). 
the corresponding exponential term will be reduced to sum of exponential of the eigenvalues of 
the operator. Similarly second exponential term can also be reduced to the sum of exponential 
of the eigenvalues of the corresponding quantum operator {dSx^^ + eSy^^ +/Sf^^). Finding 
eigenvalues of these quantum mechanical operators is not difficult; they are ±X^^^ and ±X^^^ 
respectively, where, X^^' = jVa^ + b^ + c^ and X^^^ — \\/ d'^ ^ p ■ We now write the 
transfer operator in Eq ( fT9l ) as, 

= 4 £ cosh^i^^i cosh;i^*2 (20) 

{e;-} 

where the sum runs over all the configurations of the classical spins, jO^'^'z}, within 

the quantum unit «, which is denoted by {2^}. Note, T, depends only on a^^i and o^*^- and we 
can write elements of transfer matrix P as: 

P^Fe.^Fe^ = 4 £ COSh^l^^' COSh?l^^2 g/» (21) 

where X^^^ and h are evaluated with the particular value of (o^^^ ,0^^^). Next step here 
would be to get eigenvalues of the P matrix, which in turn will give partition function for a 
particular orientation (G,(|),\|/) of the chain. Now Eq (fT6l) . Eq (fTTI) and Eq (fT5l) can be used for 
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Figure 3: yTvs.T curve for Fe2Nb. Open circles (O), connected by segmented line, represent experi- 
mental data points whereas filled triangles (A), connected by continuous line, represent calculated data 
points. The experiment was done with applied field of 100 Oe. All the theoretical points are calculated 
with same value of applied field. 

obtaining xTvs.T curve at low temperature (note, while Eq (fTTI) will give Xo^ Eq (fT?]) will give 
us Xm). 

4.1 Result and Discussion 

In Fig dll) we show the best fit xTvs.T plot with experimental data. The best fit parameters 
are: J = 20.3 cm^^ (antiferromagnetic), Nzf = 2.3x 10^^ cm^^ (antiferromagnetic), g = 2.08, 
a = 33°, (0 = 22.5° and 5 = 8°. The theoretical fit was obtained from 110 geometrical units 
or 220 chemical units. We observe that, below 6K, %r value decrease with temperature due 
to two effects: (a) finite size of chains and (b) antiferromagnetic interchain interaction. Due to 
finite size, magnetization of the chain has an upper bound and % reaches a saturation value and 
below this temperature, xT value will begin to decrease. Effect of antiferromagnetic interchain 
interaction further reduces the value of susceptibility at low temperatures. We also note that, 
at higher temperatures theoretical value of susceptibility is lower compared to the experimental 
value. Since, unlike in a real system, the direction of Fe moments in the model are fixed and can 
not change, resulting in smaller magnetization than the powder compound. Two more things 
can be estimated from our model; first is the energy for creating a domain wall in the chain, A^, 
the second being the anisotropy constant, D. 

One of the important properties of SCMs is the slow relaxation of magnetization. The 
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relaxation dynamics is controlled by the blocking temperature of the system and is related to 
the energy required to create a domain wall in the chain. The lowest energy required to create a 
domain wall is the absolute difference in energy between the ground state and an exited state, in 
which one part of the chain is magnetized in one direction and the rest of the chain is magnetized 
in the opposite direction. In the class of SCMs we are discussing, visualization of domain wall 
may be difficult due to its complex structure, especially due to the presence of both quantum 
and classical units. It is sufficient here to consider each Nb ion and the associated off-chain Fe 
ion as one quantum unit, since all of these units are identical in the absence of external magnetic 
field. In order to get ground state energy of the SCM, we solve for the eigenvalues of a quantum 
unit for different configurations of the surrounding classical units. Here the Hamiltonian to be 
considered includes interaction of Nb ion with the associated off-chain Fe ion and half of the 
interactions of that Nb ion with surrounding in-chain Fe ions. We obtain the minimum energy 
when all the three Ising spins connected to the Nb ion are in same state, say,'-i-r state. So the 
ground state of the chain is when all the Ising spins (both in- and off-chin) are in the '-i-l' state. 
Now to create a domain wall, we reverse all the Ising spins of one part, say right half, and keep 
the remaining Ising spins of left half in the initial state of '-i-l ' . This will change energy of only 
one quantum unit, where the off-chain and one in-chain Ising spins are in '-i-l' state and another 
in-chain Ising spin is in '-1' state. The difference in energy, i.e., the energy to create a domain 
wall is given by (see Fig[T]): 

= \ {Va^ + b^ + c^ - Va^ + b^ + c^ I--— i)a-.=La-'.=i 

^ 34.5K (22) 

This compares well with the experimental energy gap of 33 K obtained from log(%r) vs. l/T 
plot above temperature of 6K (to exclude effects due to finite size and interchain interactions). 

Estimation of anisotropy parameter (D) can not be done directly from the fits as our Hamil- 
tonian does not contain the D parameter. Here we employ an indirect way for its estimation. 
We have assumed that low temperature spins of Fe-ions behave like Ising spins due to large 
and negative D-value from experiment. Spin of Fe ion being 2, we have also assumed that 
spin values that can be assumed by Fe is only ±2. If we attribute the deviation of the theo- 
retical curve from the experimental curve due to an assumption of canted Ising-type spin for 
the Fe spins, then from the deviation temperature (about 12K in Fig [3]) we can estimate the 
anisotropy parameter. To do this, we solve the quantum mechanical exchange Hamiltonian, 
namely, H = JS^^ ■ S^* — D(Sf^)^, with both / and D positive, of 2-spin problem involving Fe 
and Nb spins exactly. Note, all the intrachain interactions are of this antiferromagnetic Fe — Nb 
type. We conjecture that the deviation from the experimental xTvs.T plot occurs due to sig- 
nificant population of the 3rd eigenstate relative to the second eigenstate. That is the thermal 
energy kBTj is such that the ratio of the population of the 3rd eigenstate to the 2nd eigenstate 
is l/e. This yields the result that AE = £3 — £2 ~ 3D — 27 = ksTd. Since we know / from the 
fitting, we estimate D ~ 16.3 cwT^ . This estimation is close to the experimental value of 17 
cmT^ . 
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5 Conclusion 



We have outlined a method which can be used to study thermodynamic properties of an impor- 
tant class of SCMs where alternate units are of classical and quantum nature with anisotropy 
axes of classical units being non-coUinear. We have also discussed in detail the issue of inter- 
chain interaction and showed how to incorporate it in transfer matrix technique. We applied this 
method to a real system and studied its low temperature behavior. We have also carried out an 
averaging process that allows computing the magnetic properties of a powder sample. We have 
estimated anisotropy parameter and energy associated with domain wall within this model. At 
this point it is worth mentioning that, it may be better to consider all five states (0, ±1, ±2) of 
the z-component of Fe spin instead of the two states considered here. However, in that case, 
one has to deal with little bigger transfer matrix (dimension 5x5); but the main disadvantage 
of this extension will be that it will add one more parameter to the Hamiltonian, as D will now 
appear explicitly in the model. 
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